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Q: Why would you care about a 
Python reader for the HCP MEG data? 



Bridging resources 


Dataset provider: make sure the data can be used 
with any tool. 

Software provider: make sure the software can be 
used with any dataset 

Researcher: use a baseline dataset in a most 
unconstrained way. Enabling cross-toolbox, cross¬ 
methods, cross-site comparisons 



Example 1: validation of power envelope connectivity 




Hipp et al. Nat Neuro 2012 


CTF Brainstorm HCP 
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Example 2: Benchmark against new MNE family solvers 


Tminl 


6 mne-tools.github.io/stable/auto examples/index.html#inverse-prob C 
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Compute sLORETA 
inverse solution on raw 
data 


Compute MNE-dSPM 
inverse solution on 
evoked data in volume 
source space 


Compute Rap-Music on Extracting time course 

evoked data from source_estimate 

object 


Compute source power 
using DICS beamfomer 






Extracting the time 
series of activations in a 
label 


Morph source 
estimates from one 
subject to another 
subject 


Compute cross-talk 
functions (CTFs) for 
labels for 

MNE/dSPM/sLORETA 


Compute a sparse 
inverse solution using 
the Gamma-Map 
empirical Bayesian 


Compute point-spread 
functions (PSFs) for 
MNE/dSPM/sLORETA 



Generate a functional 
label from source 
estimates 


Compute LCMV 
beamformer on evoked 
data 


Compute sparse 
inverse solution with 
mixed norm: MxNE and 
irMxNE 


Compute LCMV inverse 
solution on evoked data 
in volume source space 


Compute MNE-dSPM 
inverse solution on 
single epochs 





Challenge me! 











































































Example 3: compare machine learning tools 



• Simple and efficient tools for data mining and data analysis 

• Accessible to everybody, and reusable in various contexts 
Built on NumPy, SciPy, and matplotlib 

• Open source, commercially usable - BSD license 


Classification 


Regression 


Clustering 


Identifying to which category an object belongs 
to. 

Applications: Spam detection, Image recogni¬ 
tion. 

Algorithms: SVM, nearest neighbors, random 
forest, ... — Examples 


Predicting a continuous-valued attribute associ¬ 
ated with an object. 

Applications: Drug response, Stock prices. 

Algorithms: SVR, ridge regression, Lasso, ... 

— Examples 


Automatic grouping of similar objects into sets. 

Applications: Customer segmentation, Group¬ 
ing experiment outcomes 

Algorithms: k-Means, spectral clustering, 
mean-shift, ... - Examples 


Dimensionality reduction 

Reducing the number of random variables to 
consider. 

Applications: Visualization, Increased efficien¬ 
cy 

Algorithms: PCA, feature selection, non-nega¬ 
tive matrix factorization. - Examples 


Model selection 

Comparing, validating and choosing parame¬ 
ters and models. 

Goal: Improved accuracy via parameter tuning 

Modules: grid search, cross validation, metrics. 

— Examples 


Preprocessing 

Feature extraction and normalization. 

Application: Transforming input data such as 
text for use with machine learning algorithms. 

Modules: preprocessing, feature extraction. 

— Examples 






































Challenges we solved together 


• The ethics protocol may not allow sharing all 
aspects of the data 

• Derivative data may be stored differently from what 
your tools expects (anatomy, coordinate systems) 

• Knowing which information is which file and how to 
get the missing bits and pieces 
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Programmatic access to HCP data 
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MNE-HCP 0.1.devl2 Examples Tutorials API GitHub Page ▼ 

https: /mne-tools.github.io/mne-hcp M 


Fast interface to MEG data 


Allow us to give you a flavor by a few example queries of MEG HCP data from subject 1003007: 


# Get all entries from the MEG data header 
info = hep.read_info( 1 1003007 1 , 1 task_motor' ) 

# Get continuous MEG time series 

raw = hep.read_raw( '1003007 1 , 1 task_motor 1 ) 

# Get segmented MEG time series 

epochs = hep.read_epochs(' 1003007' , 1 task_motor' ) 

# Get all MEG time series averaged across events 
list_of_evoked = hep.read_evokeds( '1003007' # 1 task_motor' ) 

# Get details on contamination and noise sources 
annotations_dict = hep.read_annot( 1 1003007 1 , 1 task_motor' ) 

# Get precomputed independent components that compose the signal time series 
ica_mat = hep.read_ica( 1 1003007 1 , 1 task_motor 1 ) 











Run all analyses supported by MNE on HCP data 
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mne-hcp o.i.devi2 .. ./mne-hcp/auto_tutorials/index.html 



Compute forward model 



Apply reference channel 
correction 
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Run temporal 
searchlight decoding on 
event related data 



Computing ERFs from 
HCP 
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Compute inverse 
solution for evoked data 




































Load the evoked responses 
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Compute inverse solution for evoked data 




Here we’ll use our knowledge from the other examples and tutorials to compute an inverse solution and apply it on event related fields. 


# Author: Denis A. Enegemann 

# License: BSD 3 clause 

import os.path as op 
import mne 
import hep 

from hep import preprocessing as preproc 


we assume our data is inside a designated folder under $HOME 


storage_dir = op.expanduser( Wmne-hcp-data' ) 
hcp_path = op.join(storage_dir, 'HCP') 
recordings_path = op.join(storage_dir, 'hep-meg') 
subjects_dir = op.join(storage_dir, 'hcp-subjects' ) 
subject = '105923' # our test subject 

data_type = 'task_working_memory' 
run index = 0 


We’re reading the evoked data. These are the same as in Visualize evoked data 

hcp_evokeds = hep.read_evokeds(onset=' stim' , subject=subject, 

data_type=data_type, hcp_path=hcp_path) 

for evoked in hcp_evokeds: 

if not evoked.comment == 'Wrkmem_LM-TIM-face_BT-diff_M0DE-mag' : 

continue 


























Get the forward model and compute the noise covariance 
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We’ll now use a convenience function to get ourforward and source models instead of computingthem by hand. 


src_outputs = hep.anatomy.compute_forward_stack( 
subject=subject, subjects_dir=subjects_dir, 
hcp_path=hcp_path, recordings_path=recordings_path / 

# speed up computations here. Setting 'add_dist' to True may improve the 

# accuracy. 

src_params=dict(add_dist=False), 

info_from=dict(data_type=data_type, run_index=run_index)) 
fwd = src_outputs[ 1 fwd' ] 


Now we can compute the noise covariance. For this purpose we will apply the same filtering as was used for the computations of the ERF in the first 
place. See also Computing ERFs from HCP. 


raw_noise = hep.read_raw(subject=subject, hcp_path=hcp_path, 

data_type= 'noise_empty_room' ) 

raw_noise.load_data() 

# apply ref channel correction and drop ref channels 

preproc.apply_ref_correction(raw_noise) 

# Note: MNE complains on Python 2.7 

raw_noise.filter(0.50, None, method= 'iir' , 

iir_params=dict(order=4, ftype= 'butter '), n_jobs=l) 
raw_noise.filter(None, 60, method= 1 iir' , 

iir_params=dict(order=4, ftype= 'butter '), n_jobs=l) 


The default output type is "ba" in 0.13 but will change to "sos" in 0.14 
The default output type is "ba" in 0.13 but will change to "sos" in 0.14 



































Compute the inverse solution 
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iir_params=dict(order=4, 

ftype= 

'butter'), n_jobs=l) 



Out: 


The default output type is "ba" in 0.13 but will change to "sos" in 0.14 
The default output type is "ba" in 0.13 but will change to "sos" in 0.14 


Note that using the empty room noise covariance will inflate the SNR of the evkoked and renders comparisons to baseline rather uninformative. 


noise_cov = mne.compute_raw_covariance( raw_noise, method='empirical' ) 


Now we assemble the inverse operator, project the data and show the results on the fsaverage surface, the freesurfer average brain. 


inv_op = mne.minimum_norm.make_inverse_operator( 
evoked.info, fwd, noise_cov=noise_cov) 


stc = mne.minimum_norm.apply_inverse( # these data have a pretty high SNR and 
evoked, inv_op, method= 'MNE' , lambda2=l. /9.**2) # 9 is a lovely number. 

stc = stc.to_original_src( 

src_outputs[ 'src_fsaverage' ], subjects_dir=subjects_dir) 


brain = stc.plot(subject='fsaverage', subjects_dir=subjects_dir, hemi= 
brain.set_time( 145) # we take the peak seen in :ref:'tut_plot_evoked' 
brain.show_view( 1 caudal' ) # admire wide spread visual activation. 

'both' ) 
and 





























Be fascinated by visual cortex activity 



time= 145.49 ms 



















The lessons learned 


• Principle of charity: there can be good reasons why 
some data are not easily available. 

• Work together with data providers 

• Document how a data resource can be used with 
your software 

• Mind your community’s software eco-system 


Talk to each other! 



Thank you 
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